Spatiotemporal Stochastic Resonance in Fully Frustrated Josephson Ladders 



O 
O 

o 

o 

Q 



Beom Jun Kim Q , Mahn-Soo Choi h ' c , Petter Minnhagen a , Gun Sang Jeon c Q, H.J. Kim d , and M.Y. Choi d 

" Department of Theoretical Physics, Umea University, 901 87 Umea, Sweden 
h Departement Physik und Astronomie, Universitdt Basel, Klingelbergstrasse 82, CH-4056 Basel, Switzerland 

c School of Physics, Korea Institute for Advanced Study, Cheongryanri-dong 207-43, Seoul 130-012, Korea 
d Department of Physics and Center for Theoretical Physics, Seoul National University, Seoul 151-742, Korea 

We consider a Josephson-junction ladder in an external magnetic field with half flux quantum per plaquette. 
When driven by external currents, periodic in time and staggered in space, such a fully frustrated system is found 
to display spatiotemporal stochastic resonance under the influence of thermal noise. Such resonance behavior 
is investigated both numerically and analytically, which reveals significant effects of anisotropy and yields rich 
physics. 
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I. INTRODUCTION 

In a nonlinear system driven by a weak periodic signal, 
noise present in the system can operate cooperatively with the 
input signal. This effect, known as stochastic resonance (SR), 
emerges when the stochastic time scale of the noise matches 
the deterministic time scale of theJnput signal. Since its in- 
troduction two decades ago or sojd SR has been observed in 
a wide range of systems.0 Recent interest has been attracted 
to the collective spatial and temporal behavior of SR in ar- 
rays of locally or globally coupled nonlinear oscillators and in 
spatially extended continuous systems under the influence of 
noise and weak periodic driving.trQ In particular, it has been 
shown that the SR in an array of locally and linearly coupled 
bistable overdamped oscillators can significantly begnhanced 
with respect to that in a single stochastic resonator]! The op- 
timal coupling and noise strength scale with the system size, 
and the signal-to-noise ratio (SNR) in the optimal conditions 
is improved by a considerable amount. Such effects of ar- 
ray enhanced stochastic resonance (AESR) have also been 
observed experimentally in a system of coupled diode res- 
onatorsfl The AESR and other interesting aspects of the spa- 
tiotemporal stochastic resonance (STSR) are believed to have 
potential importance in the area of signal and image process- 
ing and pattern formation. 

Most of the works on STSR have been focused on model 
systems, which are relatively easy to study numerically and 
sometimes analytically (cf. Ref. For instance, widely stud- 
ied are-the one-dimensional time-dependent Ginzburg-Landau 
modeH and the Langevin-type stochastic equations for lin- 
early coupled bistable oscillators. cl On one hand, they are 
rather general in nature and applicable to various systems; on 
the other hand, relevant parameters can only be determined 
phenomenologically for a specific system. Another class of 
model worth mentioning is the kinetic Ising model.Q Although 
the Ising model per se does not have any a priori dynamics, 
it can be provided with an artificially defined temporal evo- 
lution consistent with the equilibrium properties of the Ising 
magnets. While the kinetic Ising model has been extensively 
studied mainly via Monte Carlo simulation methods in vari- 
ous contexts including STSR,Q it is not clear a priori whether 
the kinetic Ising model describes the true dynamics of a real 



physical system, especially, far from equilibrium. 

In this paper, we consider a Josephson-junction ladder 
(JJL), two one-dimensional arrays of superconducting grains 
coupled with each other via Josephson junctions, in an exter- 
nal transverse magnetic field with half flux quantum per pla- 
quette (see Fig. [I]). It has two degenerate ground states and 
switches back and forth from one ground state to the other un- 
der the influence of thermal noise and the periodic (in time) 
and staggered (in space) external current driving. Such a fully 
frustrated Josephson-junction ladder (FFJJL) has several mer- 
its as an example of the real physical system to investigate 
STSR in: The temporal evolution of the system follows, the 
well-known resistively shunted junction (RSJ) dynamics^ de- 
scribed by a set of overdamped coupled stochastic differential 
equations for phase variables on superconducting grains. Each 
parameter in the RSJ model has a clearpphysical meaning and 
can be controlled easily in experiment.!! Moreover, the JJL it- 
self has been a topic of a great number of studies due to its rich 
physic sJLrj conjunction with vortex motion and frustration ef- 
fectsjjlij Here we perform numerical simulations on the RSJ 
model, and examine the spatiotemporal behavior of SR in the 
FFJJL. Of particular interest are the effects on STSR of the 
anisotropy in the Josephson coupling along the leg- and the 
rung-direction, which turns out to yield quite rich physics. We 
also present an analytical investigation, based on the Glauber 
dynamics of vortices, and make comparison with the simula- 
tion results. 

This paper is organized as follows: In Sec. II we intro- 
duce the model and specify associated control parameters in 
the system. The model is also related to a one-dimensional 
system of vortices and to the antiferromagnetic Ising model 
in a time-dependent staggered magnetic field. A part of the 
derivation of these connections is contained in Appendix A. 
Section III presents results of numerical simulations. We first 
characterize the SR behavior of the system, and investigate 
the anisotropy effects on the STSR behavior. These simu- 
lation results are compared with the analytical ones in Sec. 
IV, obtained from the Glauber dynamics of vortices. Finally, 
Section V concludes the paper with a brief discussion about 
possible experimental setups. 
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II. MODEL 

Within the RSJ model (neglecting the capacitive and induc- 
tive effects), a Josephson junction can be characterized by the 
Josephson critical current Ij or equivalently the Josephson 
coupling J = hi j /2e and the shunt resistance R, which de- 
fine the frequency scale u>j = 2eRIj/H for the dynamics of 
the system. Throughout this paper, we always measure energy 
and temperature in units of J, time and frequency with respect 
to loj, and other induced quantities correspondingly. A JJL of 
length L in a transverse external magnetic field, depicted in 
Fig. [TJ is then described by the set of coupled stochastic non- 
linear differential equations for phase variables 



Ja sinf 



A^ + Cijit) =/,(*), (1) 

where <f>ij = <fii~4>j is the phase difference across thejunction 
(ij), N(j) denotes the set of nearest neighbors of site j, and 
Qj (t) is the random noise current on the junction with zero 
average and correlation 

(Cy(*)C<'j'(*')) = 2T [ s m,(i'j') - s (ij),(fi , )] s ( t - *')• ( 2 ) 

The Josephson coupling strength (— J^) is given by 
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a for (ij) legs, 
1 for (ij) rungs, 



(3) 



where we have introduced the anisotropy factor a be- 
tween leg.-junctions (abbreviated as legs) and rung-junctions 
(rungs). E3 The external current in Eq. (|l|) is periodic in time 
with frequency f2 and staggered in space 



Ij(t) = (-l) i+x I cos fit, 



(4) 



where each site index j has been further specified by the leg 
index £(= 1,2) and the position x along the leg, i.e., j = 
(£, x). The effects of the magnetic field are manifested in the 
gauge field Aij — —Aji, which, in the Landau gauge, takes 
the form 



A, 



for (ij) legs, 
2nfx for (ij) rungs, 



(5) 



where / = $/$o is the flux $ per plaquette in units of the 
flux quantum $o = hc/2e. 

Most of the physical properties of a Josephson junction lad- 
ders (and arrays in general) can be understood most clearly in 
terms of vortices. Accordingly, we introduce an effective vor- 
tex Hamiltonian, 



>« - / - Sf x )G(x - x')(v x , -I- 5f x ,), (6) 



which is equivalent to Eq. ([I]) within the Villain scheme (see 
Appendix A for details). In Eq. (||), the vorticity v x on pla- 
quette x is given by the directional sum of the gauge-invariant 



phase differences around the plaquette and only takes the val- 
ues ±1 while the effective interaction G(x) between vortices 
obtains the form 



G(x) = 



1 + a + VI + 2a 



(7) 



in the limit of L — > oo. The additional time-dependent flux 
$ fx (t) is due to the external driving current and given by 



Sf x (t) = 



2a - 1 
7r(l + 4a 2 ) 



(-If 7 cos m . 



(8) 



For the purpose of discussions below, one important point 
here is the following: Without external current driving (To = 
Sf x = 0) and at the maximal frustration (/ = 1/2), the 
Hamiltonian (||) has two degenerate ground state $± [see 
Fig. 1(b)] with different values of the staggered magnetization 



±- 



(9) 



The additional flux favors one of the two states ^± and thus 
lifts the degeneracy. Thus the system under consideration is 
somewhat analogous to a simple two state model and from this 
perspective the existence of a SR response may be intuitively 
expected. It should be noted here that the sign of 6f x depends 
on the anisotropy factor a, implying that the choice between 
the two states ^± depends on the value of a: For a < a° cr = 
1/2 the system favors, say, the state ^ + and for a > a® r the 
other state This can be understood from the fact that the 
external currents tend to flow through the easier links, i.e., the 
rungs for small a and the legs for large a. 

In the following sections, we will use the more complete 
and accurate model (jl]) for numerical simulations, whereas the 
effective model (^) will be the starting point of our analytic 
approach. 



III. SIMULATION RESULTS 

We first present results from the simulations of the RSJ dy- 
namics. The set of equations of motion in Eq. (|l]) is integrated 
through the use of the simple Euler algorithm with discrete 
time steps At = 0.05, and the thermal noise currents satisfy- 
ing Eq. are generated to follow the uniform distribution. 
To use the efficient tridiagonal matrix algorithm,EJ we sepa- 
rate Eq. ([!]) into two parts by change of variables (see Ref. |l2| 
for details). For the staggered input current, we set Iq = 0.1 
and O = 0.002. 

In Fig. ^, we characterize typical SR behaviors of the 
isotropic system (a = 1) with size L = 128 in three dif- 
ferent ways: (a) the time series of the staggered magnetiza- 
tion, (b) the signal-to-noise ratio (SNR), and (c) the residence- 
time distribution. In Fig. ^(a), the staggered magnetization 
m(t) = (l/L) ~J2 X (— l) x m x (t) [equivalent to Eq. (H)], where 
m x has been obtained from the plaquette sum of the gauge- 
invariant phase differences, shows clear synchronization to the 



2 



input driving current at T — 0.26. Such behavior is observed 
only in the intermediate range of temperatures. At sufficiently 
low temperatures the system does not follow the external driv- 
ing but most of the time stays in one of its two stable states 
(m(t) = ±1/2), while at sufficiently high temperatures the 
resonance effect is washed out by strong thermal fluctuations. 
The SR behavior is also clearly reflected in the power spec- 
trum defined by 



S(u) = 



dte~ lwt m{t)m(0) 



(10) 



where (• • •) denotes the average over thermal noises, and the 
values m(0) = —1/2 and 9 = 50in with tn = 2-k/Q, have 
been used in the simulations. The SNR is then defined by 



SNR = 10 log- 



in 



S 



N(Q) 



(11) 



where S is the output signal strength given by the area under 
the peak at uj = and N(£Y) is the noise level at u> = fl 
The SNR and the signal strength S [Fig. ^h)] display peaks 
at T ~ 0.26 and T ~ 0.22, respectively. It is clear from the 
time series of the staggered magnetization m(t) [Fig. §(a)] 
and also the residence-time distribution [Fig. ||(c), see below] 
that these peaks are associated with the global synchronization 
of the system to the external periodic driving. As an alterna- 
tive characterization of the SR behavior, we also measure the 
residence time L., i.e., how long the system stays in one of 
its stable states^ and show in Fig. ||(c) its distribution P(t r ) 
at various temperatures. Near T = 0.2, the distribution de- 
velops a narrow peak around t r = in/2, again signaling the 
input/output synchronization. The primary peak height P±, 
defined by the area of the distribution around t r = in/2, is 
shown in the inset of Fig. 0(c), manifesting the broad peak 
near T« 0.2. 

Near SR, synchronization of the magnetization on each pla- 
quette leads to the enhancement of the spatial correlation in 
the system, as demonstrated in Fig. 0. In Fig. |3^a), we plot the 
occupancy ratio (OR)j3 which measures how many plaquettes 
in the ladder follow (in phase) the external driving. In our case 
it is given by 



OR 



j«H-ro(t+)-ro(t_))), 



(12) 



where cos(f2trp) > and ((• • •)) represents the averages over 
time and ensemble. Note that the temperature at which the 
maximum occupancy is attained coincides with the tempera- 
ture where the SNR develops a peak. To show the enhance- 
ment of the spatial correlation in a more direct way, we plot 
in Fig. 0(b) the deviation of the average number of domain 
walls for staggered magnetization from the equilibrium value, 
An(Io) = n(0) — n(Io). Reduction in the number of do- 
main walls near the resonance temperature Tsr can be ob- 
served. Here the initial suppression of the spatial correlation 
(for T < 0.2) is due to the asynchronization of the system: 
With each plaquette out of phase by a random amount, the 



effect of the external driving is similar to that of random dis- 
order. 

We next turn to the anisotropy effects, which are summa- 
rized in Figs. |, H and H (i) The OR under goes a dramatic 
change as a is varied over a cr ~ 0.67. For a < a cr the OR 
reaches its minimum values at the resonance temperature Tsr, 
whereas maximum values are attained for a > a cr (Fig. 0). 
The OR minimum corresponds to a phase difference of ap- 
proximately ir between the applied field and the staggered 
magnetization, (ii) The SNR at Tsr is substantially sup- 
pressed near a « a cr , as shown in Fig. (iii) The resonance 
temperature Tsr in general increases with a, although it is 
ill-defined around a cr (Fig. 0; see also below). These curious 
behaviors can be understood by recalling that the anisotropy 
factor affects the system in two different ways: On one hand, 
a simulates the effective coupling strength between plaque- 
ttes [see Eq. (0)]; this is analogous to the coupling strength 
in the system of bistable oscillators, and explains straightfor- 
wardly the increase of Tsr with a. On the other hand, de- 
pending on the value of a, the system favors one of the two 
states fy± [see the discussion below Eq. (A6)], giving rise to 
the 7r-shift-like behavior in the system response (occupancy 
ratio) with respect to the signal. In addition, for a ~ a cr , the 
external currents are not efficient to drive the system since it 
is not obvious which links the currents will prefer. This natu- 
rally leads to suppression of the system response, as shown in 
Fig. [| and more clearly in Fig. Note the rather reasonable 
agreement between the crossover anisotropy a cr in simula- 
tions and a° cr = 1/2 obtained from simple consideration of 



Eq. (|A6|). 

We now concentrate on the response of the system around 
a w a cr . Surprisingly, Fig. 0(a) reveals that the SNR for 
a = 0.7, slightly above a cr , exhibits reentrance-like behavior, 
characterized by a double peak. Further, the time series of the 
staggered magnetization shown in Fig. 0(b) indicate that the 
system favors different states at lower and at higher tempera- 
tures. One possible explanation for these effects goes as fol- 
lows: For a > a cr , the currents at a given time t favor to flow 
through legs, but only slightly, with \P± almost degenerate. In 
our model described by Eqs. ([j]) and (0), the ratio of the noise 
level to the Josephson coupling on each link is given by T on 
rungs and T j a on legs, resulting in that the legs are subject to 
stronger random noise than the rungs (a cr < a < 1). As the 
temperature is raised, stronger random noises thus effectively 
hinder the currents from flowing through the legs leading to 
opposite circulation currents on neighboring plaquettes. This 
argument is consistent with the equilibrium value of the stag- 
gered magnetization (m) as a function of the temperature in 
the presence of a static staggered current (£1 = 0), displayed 
in Fig. 0(c). For those values of a, which lead to the double- 
peak structure, both temperatures corresponding to the two 
peaks have thus been plotted in Fig. H. 
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IV. ANALYTICAL RESULTS 



V. CONCLUSION 



In this section we examine the SR behavior in an analyti- 
cal way, based on the Glauber dynamicstj of vortices. Since 
the effective interaction between vortices, given by Eq. (^), 
decays exponentially with the interaction range of the order 
of the lattice constant, we can regard it essentially as the 
nearest-neighbor interaction. The effective vortex Hamilto- 
nian in Eq. (Q), together with Eq. (||), then reduces to the one- 
dimensional (ferromagnetic) Ising Hamiltonian in the oscillat- 
ing field 



Hi 



-Ji^2a x (j x+ i — h^2<7 x 



(13) 



with the (staggered) Ising spin a x = 2(— l) x m x , where the 
nearest-neighbor interaction and the oscillating field are given 
by 



Jr 



7T 2 1 + a - VI + 2a 
~2 Vl + 2a ' 



(2a-l)(l + a- yT+2^) r 

h = —2n — — ■ 1 cosllt. 

(1 + 4a 2 )Vl + 2a 



(14) 



We consider the Glauber dynamics of the Hamiltonian (Ob 
and compute the SNR of the staggered magnetization m(t) = 
(1/2L) ^2 X <j x . We follow the procedure in Ref. [l6|and obtain 
the output signal strength and the noise level 



S = 



irh 2 



^ 2 (l-7 2 ) 



ST 2 A 2 {l-~f) 2 + n 2) 



si(l - 7752) 



(15) 



where A = exp(-J z /T), 7 = tanh(2 Ji/T), 77 = 
tanh(J 7 /T), si = [(A + ifl) 2 - A 2 7 2 ] 1/2 , s 2 = A<y/(A + 
iVl + si), and SR denotes the real part. 

The SNR given by Eq. <JTT]) together with Eq. ([B]) is plot- 
ted in Fig. ||(a) for various values of a, demonstrating clearly 
the SR behavior. Furthermore, it reproduces most of the 
anisotropy effects, which have been observed in the numer- 
ical simulations. For example, Fig. |](b) reveals substantial 
suppression of the SNR at Tsr. for a w a cr , while the in- 
creasing tendency of Tsr with a is shown in Fig. ^(c). Inter- 
estingly, however, in spite of the qualitatively good agreement 
described above, one striking difference can be recognized be- 
tween the analytical results based on the Glauber dynamics 
and the simulation results obtained from the RSJ (Langevin) 
dynamics: In the former there is no reentrance-like behavior 
in the SNR and consequently no abrupt increase of Tsr near 
a w a cr . Further, the crossover anisotropy close to the value 
a° cr = 1/2 is obtained, while the simulations of the RSJ dy- 
namics have given a cr « 0.67. 



We have considered a fully frustrated Josephson-junction 
ladder, driven by external currents, periodic in time and stag- 
gered in space. In particular the spatiotemporal aspects of the 
stochastic resonance, resulting from the interplay between the 
external driving and the thermal noise, have been examined by 
means of numerical simulations performed on the resistively 
shunted dynamics. It has been found that the anisotropy in 
the Josephson coupling along the leg- and the rung-direction 
yields rich physics, e.g., resulting in a double peak SR struc- 
ture. We have also presented an analytical investigation, based 
on the Glauber dynamics of vortices, and which is in a quali- 
tative agreement with with the simulation results. 

A particularly interesting aspect of the fully frustrated 
Josephson junction ladder is the possibility of manufacturing 
such a structure and performing real measurements. Owing 
to the recent rapid advancement in lithographic technology, 
already available are large-sized Josephson-junction ladders 
of good quality. In order to verify the SR effects discussed 
in the present paper one might measure vortices on (possi- 
bly selected) plaquettes in the ladder using the superconduct- 
ing quantum interferencexlevices (SQUID) with a wide range 
of operational frequency.^ Another indirect way might be to 
measure the voltages across as many selected legs as possi- 
ble. Highly challenging from an experimental point of view 
would be to apply a staggered (but otherwise uniform) bias 
current (through each site) along the array.EJ It requires an 
impedance for each injection of current, which should be suf- 
ficiently large compared with the normal-state junction resis- 
tance. One needs to choose superconducting materials with 
higher Josephson energy, which allows higher operation tem- 
peratures and hence more choices for resistors. Another pos- 
sibility could be to use an auxiliary array of junctions, which 
provides high impedances at sufficiently high frequencies. Il3 
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APPENDIX A: 

In this Appendix, we derive the effective vortex Hamilto- 
nian (|(J). We start with the static case and consider for conve- 
nience the XY Hamiltonian: 



H = - Jy cos((f>ij - Aij), 



(Al) 
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which corresponds to Eq. in the absence of external cur- 
rentpdjOiAing (Iq — 0). Application of the duality transforma- 
tiont3 c2l at sufficiently low temperatures then decomposes the 
Hamiltonian in Eq. ( Al ) into the sum H = Hsw + Hy of the 
spin-wave part Hsw (which is given in a trivial form) and the 
vortex part 



H v = - f)G(x - x')(v x , - /). 



(A2) 



In the presence of static staggered currentsJ, ^ /: but 
SI = in Eq. (Q)], we adopt the Villain schemeEJ to obtain, to 
the linear order in Iq, the effective vortex Hamiltonian 



{v x -f- Sf x )G(x - x'){v x , - f - 5f x ,), (A3) 



where the additional effective flux Sf x depends on the fraction 
p of the external currents flowing through rungs in the zero- 
temperature configuration: 



^From Eq. (|lj) without noise current, p can be obtained 



1 



P = 



1 + 4a 2 ' 



leading to 



Sf x = 



2a - 1 
7r(l + 4a 2 ) 



(A4) 



(A5) 



(A6) 



When the external (staggered) current is periodic in time 
(SI 7^ 0), the additional flux also varies periodically in time. 
For sufficiently low-frequency driving (SI <C 1), we can take 
Eq. (^]) for the additional time-dependent flux 5f x (t). Such a 
periodic flux induces oscillations between the two states fy±, 
suggesting the possibility of the SR behavior in the system. 
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FIG. 1. (a) Schematic diagram of the fully frustrated Joseph- 
son-junction ladder (FFJJL) driven by staggered currents, (b) Two 
degenerate ground states of FFJJL in the absence of external current, 
with Q's indicating the positions of vortices. 
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FIG. 2. Various characteristics of stochastic resonance in the sys- 
tem of size L = 128 for the driving frequency Q, = 0.002: (a) 
time-series of the staggered magnetization m(t) (solid line) and the 
external current (dashed line) at T — 0.26; (b) signal-to-noise ratio 
SNR (left vertical scale) and output signal power S at frequency Q 
(right scale) as functions of the temperature; (c) residence-time dis- 
tribution at several temperatures [inset: primary peak height in the 
residence-time distribution as a function of the temperature]. 



FIG. 3. (a) Occupancy ratio OR and (b) deviation in the 
average number of domain walls from the equilibrium value, 
An ee n(0) - n(J ), for a = 1 and L =32, 64, 128, and 256. 
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FIG. 4. Occupancy ratio for different values of the anisotropy 
factor a. The lower panel (b) focuses on the region a ~ a cr and 

r ~ t sr . 
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FIG. 5. SNR, 10 log 5 (left vertical scale), and OR (right scale) at 
Tsr as functions of the anisotropy factor a. 
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FIG. 6. Resonance temperature Tsr, obtained from the SNR, as 
a function of the anisotropy factor a. It is ill-defined near a ~ a cr , 
where the system response is strongly suppressed. The two values 
of Tsr for a > a cr correspond to the double-peak structure of the 
SNR (see Fig. m. Behavior of Tsr from the occupancy ratio as well 
as from the output signal strength is shown in the inset. 
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FIG. 7. Reentrance-like behavior near ot cr : (a) the SNR as a 
function of the temperature for a — 0.7 and (b) the corresponding 
time-series of the staggered magnetization m(t). For clarity, vertical 
positions of the data curves are shifted with respect to each other. 
In (c), also shown is the equilibrium value of the staggered mag- 
netization as a function of the temperature in the presence of static 
staggered currents (fl = 0). 
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FIG. 8. Resonance behavior displayed by the Glauber dynamics 
of vortices in FFJJL: (a) SNR as functions of temperature for various 
values of a; (b) SNR at Tsr as a function of the anisotropy factor a; 
(c) Tsr as a function of the anisotropy factor a. 
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